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Abstract 

Hard spheres with an attraction of range a tenth to a hundredth of the sphere diameter are constrained to 
remain fluid even at densities when monodisperse particles at equilibrium would have crystallised, in order 
to compare with experimental systems which remain fluid. They are found to have a fluid-fluid transition 
at high density. As the range of the attraction tends to zero, the density at the critical point tends towards 
the random-close-packing density of hard spheres. 

PACS: 82.70.Dd Colloids, 61.20.Gy Theory and models of liquid structure 



Argon forms a liquid because argon atoms attract 
each other and these dispersion attractions between 
the atoms are relatively long-ranged; the volume over 
which one argon atom attracts another is comparable 
to the volume one argon atom of the pair excludes to 
another. If we could reduce the range of the attrac- 
tion between argon atoms then the liquid phase would 
disappear from the equilibrium phase diagram when 
the volume over which the atoms attract was of order 
one tenth of the volume they exclude to each other. 
Of course we cannot change the interaction between 
argon atoms but there are well-established colloidal 
systems whose interactions we can change. The liq- 
uid phase disappears from the equilibrium phase di- 
agram because the fluid-fluid transition is preempted 
by the crystallisation of the fluid. But although the 
fluid-fluid transition has disappeared from the equilib- 
rium phase diagram of monodisperse particles, experi- 
ments often do not observe crystallisation, presumably 
due to a combination of a large free energy barrier 
to crystallisation and the destabilising effect of small 
amounts of polydispersity on the crystalline phase. 
As crystallisation does not occur it does not preempt 
the fluid- fluid transition, which is therefore observ- 
able. With this in mind we study the behaviour of 
spherical particles with a short-range attraction which 
are constrained to remain fluid. We study attraction 
ranges down to a hundredth of the diameter of the 



hard core — this is what we mean by very short-range 
attractions. We find that as the range decreases, the 
density at the critical point increases to very high val- 
ues. For a sufficiently short range the critical point 
lies above the density of the kinetic glass transition 
observed in experiments on hard-sphere-like colloids. 
Experiments on colloids with very short-range attrac- 
tions have found a threshold beyond which diffusion 
ceases 0,0, and Poon, Pirie and Pusey || have pre- 
viously suggested that this is due to an arrested fluid- 
fluid phase separation. If as we find, the dense fluid 
has a density above that of the hard-sphere glass tran- 
sition, then it is not surprising that the dynamics of 
phase separation become arrested. 

Here we will not consider the crystalline phase at 
all. Our results are for a system of particles which 
is constrained to remain fluid at all temperatures and 
pressures; see refs. [|]-|| for a discussion of the appli- 
cation of constraints to stabilise a phase which would 
otherwise be metastable or unstable. Although exper- 
iments on near-monodisperse colloidal spheres show 
that they crystallise readily, at least as long as the 
attraction is not too strong, polydisperse colloidal 
spheres often never crystallise |?J and the presence 
of a very short-range attraction makes the crystalline 
phase even more sensitive to polydispersity ||. By 
polydisperse spheres we mean that the spherical par- 
ticles do not all have the same diameter but have a 



range of diameters. Our theory is a perturbation the- 
ory about a hard-sphere fluid and so completely ne- 
glects the crystal. Thus we will not need to explicitly 
apply a constraint within the theory. We do, however, 
need to assume that it is possible to apply a constraint 
to the system which has almost no effect on the fluid 
phase but completely prevents crystallisation. 

We chose a simple potential with a hard-sphere 
core and an attraction in the form of a Yukawa func- 
tion. The hard-sphere-l- Yukawa potential is a spher- 
ically symmetric pair potential so the interaction en- 
ergy v depends only on the separation r of the centres 
of the two particles, 



v(r) 



DC 



-e£exp[/c(l -r/a)} 



r < a 
a < r 



(1) 



where a is the hard-sphere diameter and e is the en- 
ergy of interaction for touching spheres. With this 
potential the thermodynamic functions depend on the 
reduced temperature kT/e, and the reduced density 
7] = (N/V)(7r/6)a' i which is the fraction of the vol- 
ume occupied by the cores of the particles, k, T, N 
and V are Boltzmann's constant, the temperature, the 
number of particles and the volume, respectively. 

We require a free energy for this potential which is 
accurate up to very high densities, up to near random 
close packing which is at a volume fraction rj ~ 0.64- 
0.65 |||nj. Speedy [|10| has obtained, from computer 
simulation data, an accurate equation of state of hard 
spheres up to random close packing. This enables us 
to use a perturbation theory, i.e., to start from the 
Helmholtz free energy in the infinite temperature limit 
of our model, which is hard spheres, and add on the 
energy as a perturbation. Then our expression for the 
Helmholtz free energy per particle a at a temperature 
T and a volume fraction 77 is 



pa(v,T)=0a hs (v)+Mri,T), 



(2) 



where a^ s is the Helmholtz free energy of hard spheres, 
u is the energy per particle and (3 = 1/kT. As the en- 
ergy of a fluid of hard spheres is zero, j3a,h s = —Sh s /k, 
where Sh s is the entropy per particle of hard spheres, 
which is, according to Speedy flC, 11] 



f|£ = l-lnp+CHrio-ri + So+N-HnNgim), (3) 



where C = 2.8, So = 
N g(Vo) = exp 



-0.25 and N g is 



N (a- 7(7/0 - Vm)' 



(4) 



where a = 2, 7 = 193 and r\ m = 0.555. In these 
equations the value of 770 at any density is determined 



by minimising the free energy at that density. This 
form of the free energy is optimised for the dense fluid. 
Essentially, if we start from any configuration of the 
dense fluid and begin to expand all the spheres (so 
increasing the volume fraction) then at some point 
the spheres will touch and then the spheres cannot be 
expanded further. At this point the volume fraction 
is 770; this can be seen from the log term in eq. (||) 
which diverges when 77 = 770- If we start from differ- 
ent configurations then after expansion of the spheres 
we may end up with a different value of 770 • The 
larger the difference 770 — r\ then the more room the 
spheres have which increases the entropy. However, 
simulations show that there few arrangements of the 
spheres which have a large 770, therefore there is an 
entropic cost to being in an arrangement with a large 
770. N g (r]o), eq. (||), is essentially the number of ways 
of arranging spheres such that the maximum possible 
volume fraction is 770; it is maximal at 770 = r] m . The 
competition between the third and fifth terms in eq. 
(||) then determines the value of 770. As the spheres 
touch when 77 = 770 and if we assume that the expan- 
sion is isotropic then the separation b of spheres at a 
given 77 and 770 is 



b/cr = (voh) 1/3 , 



(5) 



just as in a crystal. The energy of attraction is ap- 
proximated by the energy of interaction of each sphere 
with its six neighbours § at a separation b 



u = -3v(6) = -3t;((7/o/7/) 1/3 ). 



(6) 



As the energy depends on 7/0, the total free energy, 
eq. (||), is minimised to obtain 7/0 at each density and 
temperature. 

Guides to the accuracy of our free energy are ob- 
tained by comparison with existing simulation data. 
For k = 7, Hagen and Frenkel |12| find a fluid-fluid 
critical point at kT/e = 0.41, 77 = 0.26, whereas we 
predict kT/e = 0.54, 77 = 0.30. The agreement is fair 
although not quantitative and we expect our theory 
to do better at higher densities. Applying an approx- 
imation of the type, eq. (||), to a face-centred-cubic 
(fee) crystal |13|] yields an fcc-crystal-fcc-crystal crit- 
ical point when k = 100 at kT/e = 1.1, 77 = 0.69. 



Bolhuis, Hagen and Frenkel 14] using computer sim- 
ulation and perturbation theory predict kT/e = 0.70, 
77 = 0.71. Again there is fair but not quantitative 
agreement 



Results for four, short, ranges are plotted in fig. [j] 
[11]. A simple liquid such as argon is reasonably well 
modeled by an attraction of inverse range n = 1.8 fig] . 
The results are for inverse ranges up to two orders of 
magnitude greater. The notable feature is that the 
critical densities and the densities of the liquid phase 
are high and move to higher density as the range de- 
creases. At high density the particles are pushed to- 
gether until they are within range of the attraction. 
This occurs at separations between the surfaces of the 
spheres b — a = 0(aK^ 1 ). With the particles just 
within range of the attraction there is a clear ener- 
getic driving force towards phase separation: the fluid 
lowers its energy at fixed overall density by some of 
the fluid condensing into a dense fluid where all the 
spheres are well within the range of the attraction of 
their nearest neighbours. This is just what was ob- 
served by Bolhuis, Hagen and Frenkel [O] in the fee 
crystal. In the absence of the crystalline phase, due 
to polydispersity perhaps, the transition simply shifts 
over to a fluid-fluid transition and it occurs at a lower 
density due to the fact that the random close-packing 
density which is the maximum density of amorphous 
spheres is lower than the maximum density of spheres 
achievable in an fee crystal. Because of the smaller 
number of neighbours in the dense fluid as compared 
to the crystal the transition shifts to a lower tempera- 
ture but in both cases the critical temperature varies 
little with changing range. 

Grant and Russel |l]], and Verduin and Dhont 
ly have performed experiments on colloids which are 
hard-sphere-like at high temperature but as the tem- 
perature is reduced the solvent becomes a poor solvent 
for the alkane chains which are grafted to the surface 
of the colloids. There is then a very short-range at- 
traction when the grafted layers of two colloids over- 
lap. Verduin and Dhont estimate the range of the 
attraction to be less than lnm for colloids with a di- 
ameter 80nm. They assess the strength of the attrac- 
tion via a parameter r [ITJ, 17 1 which is related to the 
second virial coefficient B 2 by 



- (l - B 2 /B 



(7) 



where B 2 is the second virial coefficient of hard 
spheres. We have replotted the phase diagrams for 
k = 20 and 100 in the density-r plane in fig. [2|. The 
experimental results are in the range r = 0.1 to 0.2, 
and r\ = 0.1 to 77 = 0.4. These densities and temper- 
atures lie within the coexistence region for k = 100. 
They locate both a spinodal and a 'static percolation' 
line where diffusion of the particles stops. At these 



densities and temperatures, the fluid is unstable with 
respect to phase separation into a more dilute phase 
and a very dense fluid phase — this phase has a vol- 
ume fraction ~ 0.56. Beyond a volume fraction of 
0.56-0.58 1 18], the relaxation time of a fluid of hard 
spheres exceeds typical experimental times which are 
of the order of 1000s. This is referred to as a ki- 
netic glass transition, above it the samples are not 
equilibrated on the experimental time scale and so 
are not an equilibrium fluid. Thus as phase separa- 
tion proceeds and domains of this dense phase appear 
the phase separation dynamics may become arrested 
due to the very slow relaxation within these dense do- 
mains. This is the scenario suggested by Poon, Pirie 
and Pusey on the basis of their experiments on colloids 
with a longer range but still short-range attraction ||] . 
Alternatively, as Grant and Russel []l| have suggested 
the phase separation may be fluid-crystal phase sepa- 
ration. 

We predict a critical point at a density which in- 
creases as the range decreases and at a value of r 
which also increases as the range decreases. Although 
we cannot perform calculations at zero range, k = 00, 
extrapolation of our results together with the results 



of Bolhuis, Hagen and Frenkel [14] who were able to 
study the zero range limit in the crystal, suggests that 
in the zero-range limit there is a fluid-fluid critical 
point at the random-close-packing density. The criti- 
cal point would be at roughly the same temperature 
as the critical points in fig. [I], which implies that 
as in the crystal [14| it occurs at infinite r. Baxter 
Jl~6| ] solved the Percus-Yevick (PY) approximation for 
hard-spheres with a zero-range attraction. Within the 
PY approximation there are 2 routes to the thermody- 
namic functions. If the compressibility route is used 
the critical point is at the low volume fraction 77 = 0.12 
and at r = 0.098 []19| , whereas via the energy route the 
prediction is 77 = 0.32 and r = 0.12 [19|. Our results 
suggest that the critical point predicted by the PY ap- 
proximation may be an artifact of this approximation. 
Stell [17] has shown that the virial expansion is patho- 
logical for all finite r in the limit of the range of the 
attraction tending to zero. However, this pathology 
originates in crystalline clusters which we have elim- 
inated with our constraint that crystalline configura- 
tions are not allowed Q || . So the pathological virial 
expansion shows that in the zero-range limit, the equi- 
librium, unconstrained, fluid is completely unstable at 
finite r; see Refs. [14, 1?], |2[J, |2l|] for the equlibrium 
phase diagram in the zero-range limit. 

In conclusion, we have determined the phase dia- 
gram of hard spheres with an attraction with a range 



of order 0.1 or 0.01 of the hard-core diameter, which 
are constrained not to crystallise. The fluid- fluid tran- 
sition persists, according to our approximate theory, 
for all ranges of the attraction. As the range decreases 
the density at the critical point increases and can be- 
come very high, near the random-close-packing den- 
sity of hard spheres. As the density is so high observ- 
ing it will be difficult as the dynamics are very slow 
at these densities; the densities can exceed that of the 
glass transition of hard spheres. Due to these slow 
dynamics a glass-glass transition may be observed in- 
stead of a fluid-fluid transition but again it may be 
impossible to observe directly. The difficulty in ob- 
serving fully equilibrated coexistence does not mean 
that the transition has no observable consequences. 
Out of equilibrium systems tend to head toward equi- 
librium and even if they do not reach equilibrium their 
final state may be, roughly speaking, the point on the 
path to equilibrium where the dynamics stop. This is 
our tentative interpretation of the results on colloids 
with a very short-range attraction |l|,|2]]: the ceasing 
of diffusion observed is arrested fluid-fluid phase sep- 
aration. One final point is that as the range decreases 
the critical point, with its associated large fluctuations 
and critical slowing down of the dynamics |^2[ , will ap- 
proach the kinetic glass transition. What effect this 
will have on the kinetic glass transition is unknown. 
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Figure 1: The fluid phase diagrams in the temperature-density plane, for four different ranges. The curves 
are, from left to right, for inverse ranges k = 7, 20, 40 and 100. In each case the curve encloses the fluid-fluid 
coexistence region and the critical point is the highest point on the curve. 
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Figure 2: The fluid phase diagrams in the r-density plane, for two different ranges. 




